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Abstract 

We derive an asymptotic solution of the Einstein field equations 
which describes the propagation of a thin, large amplitude gravitational 
wave into a curved space-time. The resulting equations have the same 
form as the colliding plane wave equations without one of the usual 
constraint equations. 

1 Introduction 

Gravitational waves are one of the most important features of Einstein's 
general theory of relativity. The Einstein field equations are highly nonlin- 
ear, and a question of fundamental interest is how nonlinearity affects the 
propagation of gravitational waves. 

Small amplitude gravitational waves are well described by the linearized 
Einstein equations which completely neglect nonlinear effects. Large am- 
plitude unidirectional gravitational plane waves are described by the exact 
Brinkmann- Rosen solution of the vacuum Einstein equations ||], |2| . Despite 
the nonlinearity of the Einstein equations, a gravitational plane wave prop- 
agates into flat space-time without distortion, and there are no dynamic 
nonlinear effects. 

The simplest situation in which nonlinear effects are significant is when 
a large amplitude gravitational wave propagates into a curved space-time. 
An important special case is when the space-time ahead of the wave is that 



of a counter-propagating gravitational plane wave. The resulting space-time 
has a two-parameter family of spacelike isometries, and the metric is given 
by the exact colliding plane wave solution of the vacuum Einstein equations 
H |], H ||]. In the case of more general space-times ahead of the wave, exact 
solutions do not exist. 

In this paper, we derive an asymptotic solution of the Einstein equa- 
tions which describes the propagation of a thin, large-amplitude, pulse-like 
gravitational wave into a general curved space-time. The solution applies 
provided that the metric varies much more rapidly inside the wave than 
on either side of the wave. As a result, the wave can be approximated 
locally by a nonlinear plane wave which is distorted as it propagates into 
the curved space-time. For plane-polarized waves, the asymptotic solution 
is given by equations flOl) , (^6|) , (|3?7|) , and (|3^)-(^Tl|) below. For non- 
polarized waves, the asymptotic solution is given by (|3.1| ), ( |6.3j ), (|6.5D, and 
(|6.8[) -( |Ol|) . The asymptotic equations consist of the colliding plane wave 
equations without one of the usual constraint equations. The colliding plane 
wave equations are therefore canonical equations for nonlinear gravitational 
waves, and they describe a much larger class of solutions than the ones with 
exact plane-wave symmetry. 

The nonlinearity of the asymptotic equations may lead to the develop- 
ment of a space-time singularity. A plane gravitational wave propagating 
into flat space-time does not steepen. Consequently, the mechanism of singu- 
larity formation in gravitational waves differs from the nonlinear steepening 
of waves in quasilinear hyperbolic systems which leads to the formation of 
shocks. Instead, the singularity formation is caused by the mutual focusing 
of the gravitational wave and the curved space-time into which it propa- 
gates. A second nonlinear effect described by the asymptotic equations is 
the permanent distortion of space-time by the passage of a gravitational 
wave. A curved pulse generates a backscattered gravitational wave which 
propagates into the space-time behind it. 

In Section 2, we briefly summarize the exact colliding plane wave solution 
of the Einstein equations. In Section 3, we give an overview of the asymptotic 
expansion. In Section 4, we write out expansions of the metric components, 
the connection coefficients, and the Ricci curvature components. In Section 
5, we construct a coordinate system in which the metric adopts its simplest 
form. In Section 6, we complete the derivation of the asymptotic equations. 
In Section 7, we show that the same equations follow from an expansion of 
the variational principle for the Einstein equations. In Section 8, we explain 
how to derive boundary conditions for the asymptotic equations, and in 
Section 9 we consider some specific physical examples. 
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2 Colliding plane waves 

The vacuum Einstein field equations imply that 

Ricci = 0, (2.1) 

where Ricci is the Ricci tensor associated with the metric tensor g. The 
plane-plolarized, colliding plane wave solution of fl2.ip is given by 

g = -2e~ M du dv + e~ u {e v dy 2 + e~ v dz 2 ), (2.2) 

where the functions M(u, v), U (u, v), V(u, v) satisfy the colliding plane wave 
equations, 

(2.3) 
(2.4) 

(2.5) 

(2.6) 

(2.7) 

Equations (|2.3|)- (l2.5| ) are wave equations for M, U and V in characteristic 
coordinates (u, v). Equations (^)-(|2?7|) are constraints which are preserved 
by (fO|)-(2.5). To specify a unique solution, the wave equations can be 
supplemented by characteristic initial data for M, U, V on the lines u = 
and v = which satisfy the appropriate constraint equations. 

The metric which describes the collision of non-polarized plane waves is 

g = -2e' M du dv + e~ u (e v cosh Wdy 2 - 2 sinh Wdy dz + e~ v cosh Wdz 2 ), 

where the functions M(u, v), U(u,v), V(u,v), W(u, v) satisfy 



U U v — 


U U U V , 


'uv — 


~ (U U V V + U V V U ) , 


M uv = 


i (-u u u v + v u v v ) , 


U uu — 


\{U 2 U + V 2 )-U U M U , 


U vv — 


\{U 2 + V 2 )-U V M V . 





— U u U v , 


(2.8) 




= ~(U U V V + U V V U ) - (V U W V + V V W U ) tanh W, 


(2.9) 


w uv 


= \ (U U W V + U V W U ) + V U V V sinh W cosh W, 


(2.10) 


M uv 


= ^(-U U U V + V U V V cosh 2 W + W U W V ), 


(2.11) 




= \{U 2 U + V 2 cosh 2 W + W 2 )- U U M U , 


(2.12) 




= \{U 2 + V 2 cosYv 2 W + W 2 )-U v M v . 


(2.13) 
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When W = 0, this solution reduces to the plane-polarized solution. When 
all functions are independent of v, the solution reduces to the Rosen form 
of the exact unidirectional plane wave solution. 



3 Overview of the expansion 

In this section, we outline the main ideas of the derivation of the asymptotic 
solution. For simplicity, we describe the case of plane-polarized waves. The 
algebraic details are given in the following sections. 
We consider metrics of the form 

/ u(x) \ 

g = g(-^Wj, (3.1) 

g(9,x;e) = g(9,x)+ek(9,x)+0(e 2 ), 

where e is a small parameter and u is a scalar-valued phase function with 
du ^ 0. This ansatz corresponds to a metric that varies rapidly and strongly 
in the u-direction. The phase u is a null function of the metric, at least up 
to the order e. That is, it satisfies 

g 9 {du,du) = 0(e 2 ), (3.2) 

where g" is the contravariant form of the metric tensor. The component 
form of this equation is written out in Q4,5| ) below. 
The scaled variable 

e = H (3.3) 

is a "stretched" coordinate inside the wave. We assume that the derivatives 
of g(9,x;e) with respect to 9 decay to zero sufficiently quickly as 9 — > oo. 



Thus, the solution (3J.) represents a thin, pulse- like gravitational wave lo- 
cated near the null surface u = 0. For example, if the metric is independent 
of 9 when \6\ is sufficiently large, then the solution represents a thin "sand- 
wich" wave which separates slowly varying metrics on either side. 

The Ricci tensor associated with the metric ( |3.1| ) has an expansion of 
the form 

1-2 i-i 
Ricci = — Ricci +- Ricci +0(1). (3.4) 

At leading order in e, the Einstein equations ( [2.1| ) imply that 

-2 

Ricci= 0. 
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This equation is a nonlinear, second order ordinary differential equation in 
dg for the leading order term of the metric in which the "slow" variables x 
occur as parameters. We write it symbolically as 



N(dj 



0. 



(3.5) 



In suitable coordinates (u,v,y,z), a solution of this equation is the plane- 
polarized plane wave metric 



o 



2e~ M dudv + e~ u (e v dy 2 + e~ v dz 2 ) 



(3.6) 



where M, U, V are functions of (9,v,y,z). For a metric of the form (j^ 
equation (|3~1]|) reduces to the ^-constraint equation, 



[U^ + Vi) -U e M e . 



(3.7) 



At the next order in e, the Einstein equations imply that 

-i 

Ricci= 0. 



This is a linear equation for g of the form 



F(dg,d v ,d y ,d z 



(3. 



where L is a second order linear ordinary differential operator in 8q acting on 
1 o 

g, with coefficients depending on g, and F is a nonlinear partial differential 
o 



operator acting on g. The equations in (3.8) are not independent. The 

requirement that ( jjjlf ) can be solved for g implies that M, U, and V satisfy 
the equations 



Vev = \ {U e V v + U v V e ) , 

Mov = o (-UeUv + V e V v 
2 



(3.9) 
(3.10) 

(3.11) 



Equations (pTEj )— ( 3.11 ) are identical to the evolution equations (|2.8|) — ( 2.10 ) 
for the exact colliding plane wave solution, with 6 = u/e. The leading order 
solution satisfies the constraint equation ( |3.7|) in the "fast" phase variable 9, 
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but it need not satisfy the constraint equation ( |2.7| ) in the "slow" variable v. 
If the v-constraint equation does not hold, then the asymptotic expansion 
of the metric contains higher order terms which are absent in the exact 
colliding plane wave solution. 

Equation ( |3.6[ ) implies that d v = — e M g" • du. Thus, d v is a vector on 
the light cone which is tangent to the null surface u = 0, and the "slow" 
derivative with respect to v which appears in (|3.9| )-( j3Tl"|) is a derivative 
along the bicharacteristic null geodesies associated with u. The transverse 
variables y and z occur as parameters. Therefore, in the short-wave limit 
considered here, the (1 + 3)-dimensional field equations reduce to (1 + 1)- 
dimensional asymptotic equations along the set of null geodesies associated 
with the phase u. The parametric dependence of the solution on y and 
z allows the pulse to be compactly supported in the transverse directions, 
so that the wave need not have infinite extent. Moreover, the asymptotic 
solution need not have any special exact symmetries. 

The asymptotic equations for non-polarized gravitational waves are ob- 
tained in a similar way. They consist of the general colliding plane wave 
equations fl2,g|)-( 2~T2^ ) with u replaced by 9. The v-constraint equation ( |2.13| ) 



is not required to hold. 

Since the asymptotic equations follow from the order e~ 2 and order e _1 
components of the field equations, the asymptotic solution remains valid in 
the presence of matter with a slowly varying, order one energy-momentum 
tensor, T = T(x). 

One subtle point in carrying out the expansion concerns the choice of 
the phase function u. In order for ( |3.5| ) to have a nontrivial solution, the 
phase u must be a null function of the leading order metric, but u need not 
be a null function of the entire metric. However, it follows from the analysis 
in Section [| that we can use a transformation of the form 

m£*^,j;;e) (3.12) 

to choose a phase which satisfies ( |3.2[ ), The asymptotic solutions obtained 
with the use of the old and the new phases can be shown to be equivalent. 
When the phase satisfies ([T^), variations in the metric propagate along 
the null geodesies associated with the phase, and the asymptotic equations 
adopt their simplest form. 
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4 Expansion of the metric and the curvature 

In this section, we write out expansions of the metric components, the con- 
nection coefficients, and the Ricci curvature components. 
We use local coordinates x a in which 

g = g af3 dx a dxK (4.1) 

Here and below, Greek indices a, (3, fj,, u, . . . take on the values 0, 1, 2, 3. We 
look for an expansion of the metric components as e — > of the form 

( U ( X ) \ fAC\ 

9a(5 = 9a/3 ,X]£ , (4.2) 



g a/3 (9,x;e) = 9 a/3 (9, x) + e g a/3 (9, x) + 0(e 2 ). 
The contravariant metric components g a @ satisfy 

Expansion of this equation in a power series in e gives 

g af3 =9 aP -e9 aP + 0(e 2 ). (4.3) 

■ — . o „ o 
In (|4.3|), 9 ap is the inverse of 9 a p, and we use the leading order metric 

components to raise indices, so that 

g^Jg^gP- . (4.4) 

With this notation, the order e term in the expansion of the contravariant 

metric component g a @ is — 9 a ^ , not g a @ . 

In terms of the metric components, we have 

g»(du,du) = ff>*L*L 



dx a dx@ 

-D du du i a du du 2 . / A . 

9 7: — ^ — — r + O e 2 . (4.5 

dx a dxP dx a dxP ; y ' 



Thus, the null condition ( |3.2[ ) holds provided that 



„ fl du du 1 nR du du , , , 

9 ^ = Q p = (4 6) 

dx a dx? ' dx a dx? K ' 
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The first condition in (4.6) states that u is a null function of g. The second 
condition implies that the phase is a null function of the perturbed metric, 
at least up to the first order in e. 

The covariant components R a g of the Ricci tensor are given by 



Roc/3 



dx x 



dx c 



where T x a p are the connection coefficients 



1 „a m ( d 9/3^ dgap 



dg 



a/3 



(4.7) 



(4.8) 



^From 
given by 



where 



dx a dx@ dxV- 
the derivative of a function f a p(,d,x), with respect to x^ is 

dfa/3 1 



dx^ 



;faf3,0Ufj, + fa0,(i, 



(4.9) 



du 

dx^ ' 



f 



a/3,0 



df 



a/3 



oe 



9 fat 



a/3,(i 



dxt 1 



We use (g^), (0), and |L9|) in (g7|) and Q) and expand the result with 
respect to e. After some algebra, we find that 



A 1 _1 a A 

r a/3 = ~ T a /3 + T a/3 + 0(e), 



ft 



a/3 



1-2 1-1 
"2 fta/3 +7 fta/3 +0(1), 



£ 



(4.10) 



where 



1 a/3 
°A 

1 a/3 



1 

2 9 

1 o 

^ 9 





S73 M ,0 u a + g a ^ e up- g al3>e Up I , 



ooo 

S/3/i,a + 9 a fj,,/3 ~ 9 a /3,(L 



1 

+ 2 9 

1 i 
-- 5 



X(x 



X/x 



111 

#/3/Lt,0 U a + g afli e u (3— 9 a f3,e Ufx 



9/3(1,9 U a + g a ^e Up- 9 a f3fi , 



(4.11) 



-2 
ft a/3 
-1 

ft a/3 



-1 



-1 -1 -1 -1 -1 

r ^ a /3,eU(j,— r M i3(tfiu a + r %/3 r "/«/— r r 
-10 

/'...., 



-1 



r ^a/3,(i— T /i /3/i,a+ t^af3fiU( 1 — f A * j3(i,9Uc 

-1 -1 -1 -1 

, -p /i -p^ _j_ -p£i -p f _ r (l r V _ -p£i -p V 

T 1 a/3 1 (iv \ 1 a/3 1 ^ti/ 1 of 1 (3(i 1 a;/ 1 /3/x- 



s 



The component form of the field equations (^^) is 

Rap = 0. (4.12) 



Using ( |4.10| ) in ( 4.12j ) and equating coefficients of e 2 and e 1 to zero, we 
get that 

Ra/3 = 0, (4.13) 

Ra/3 = 0. (4.14) 

In order to solve these equations, we first use a coordinate transformation 
to simplify the form of the metric. 

5 Coordinate transformations 

In this section, we show that there is a choice of a local coordinate system 

x a in which u = x° and the metric has the form 

o n , o , 
g = 2 001 dx u dx 1 + 9 ab dx a dx b 

+e 1 2 g la dx 1 dx a + 9 ab dx a dx b ^ + 0(e 2 ). (5.1) 

Here and below, indices a,b,c, . . . take on the values 2, 3, while indices 
k, . . . take on the values 1, 2, 3. 

The corresponding expansion of the contravariant form of the metric 
tensor is 

g S = 2g 01 d d 1 +g ab d a d b 

-e (2 g m g ab g lb d Q d a + g ac g bd L d a d b ) + 0(e 2 ). (5.2) 



For this metric, we have 



> = , g°° = o. (5.3) 



Thus, the phase u = x° satisfies (4.6), and hence (|3.2|). 

The most general coordinate transformation which is compatible with 
the expansion Q4.2]) has the form 



— -> *° ( —,x J +e q> °f —,x j +0(e 2 ), (5.4) 
x i ^( x )+ei l (^,x\ +e 2 l l (^,x) +0(e 3 ). (5.5) 
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We suppose that the phase is given by u = x° in both the old and the new 
coordinates. Thus, the change of coordinates fl5.4| ) implies a change in the 
phase of the form (3.12). 

First we simplify the leading order metric components by means of a 
transformation of the form 



x u -> x u , x l ^x % +e^ l [ — ,x . (5.6) 



x 







Expansion of the transformation law for the change in covariant tensor com- 
ponents implies that the leading order metric components transform under 
(|5.6|) according to 

1,0 1,1,0 

9oo -» 9oo+2^ k g g k + ^;e^,e9kh (5.7) 

1,0 

%i Soi + tf^Sw, (5.8) 
% -» % . (5.9) 

If the matrix 9%j is non-singular, then fl5,8| ) implies that we can transform 

9oi to zero. This contradicts the requirement that a; is null (cf. 0, Section 
109). Hence, we must have 

det 9ij= 0. (5.10) 

By an appropriate renumbering of the i-coordinates, we can suppose without 
loss of generality that 

det 0. (5.11) 

^From (|5.7j ) and ( |5.8[) , we can then choose the transformation (|5.6|) so 
that 

ooo 

9oo=9o2 =9 03 = 0. (5.12) 

. _. o 

Solving equation (|5.10|) for 5 11, we get 

9n=9 ab 9 la 9i b , (5.13) 

o , o 
where 9 is the inverse of 9 a b- We define 



9 ab g u . (5.14) 
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iFrom (glg)-(glg), it follows that 

c d 
9u=9cd g 9 , 







g\ a =g ac g ■ 



(5.15) 



Using (5.12)-(5.15) in ( |4.1| ), we find that in the transformed coordinate 
system, the metric has the form 



g = 2 Soi dx°dx 1 + g a b {dx a + g a dx 1 ) [dx b + g b dx 1 ^j + 0{e). 
^From (|4.13|), the metric (|5.16|) must satisfy the condition 



(5.16) 



Rab = 0. 



(5.17) 



Using ( [Up in (p~Tl]) , we find that 

-2 1 o 



B„<- -^{9 m ) 2 g a c g c fi %d g% 



(5.18) 



Equations Q5.17D - Q5.18D imply that 

= o, 

so g a is independent of 9. This fact allows us to remove g a by a transforma- 
tion 

s -^* ^ 1 ,^). (5.19) 



The form of the metric ( 5.16 ) is unchanged by ( 5.19 ), and 



o 

g a b 



-1 \a I „c 



(A-r c (g c + ^i), 



(5.20) 



where (A%) = (^ c ). ^From ( |5.20| ), we can set g a = 0. The metric flO 
then reduces to 



g = 2 %i dx°dx 1 + g ab dx a dx b + 0(e). 



(5.21) 



Next, we simplify the form of g. We consider the transformation of 
coordinates 



X — > £ $ 



I n / X' 



-,a; 1 +e 2 *° 



x i + e 2 P (— ,x ) . (5.22) 
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Under the action of ( 5,22; ), the form ( 5.2 1] ) of the metric is unchanged at 
order zero and 



o i o 
Soi-^je Soi • 

At order one, the components transform according to 

1 1 n 1 1 n 2 n 

001 -> *,fl^01+(*,o+*,fl)5oi, 
1 1 n 1 2,0 

00a -> *,e5oa + *?0 0a6, 
1 1 10 

0n -> 0ii+2^^ffoi, 

1 1 1 „ 

01a -> 01a + *, o 0Ol, 

1 1 
0afe — ► 0a6 ■ 



These transformations can be used to make 

0. 

The resulting metric then has the form given in (|5 



011=0Oa = 



(5.23) 



Use of fl5,l| ) and ( |5.2[ ) in (4.11) implies that the nonzero connection 
coefficients at the orders e~ l and e° are 



1 oo —0 0Ol,i 



r 1 

1 ab 



1 o 01 

2 S 0a6,l 



-l 1 o ac 

r a 0b = 0fec,0 ; 



1 00 —0 001,0) 1 0a — 7^ (0Ol,a + Sla,e) ~ 7j ^ ^abfii 



1 ab 
°1 

r afe 



t-i a 
1 06 





r i6 



1 o 01 




2 

1 o 01 1 

~2 5 (5a6,O+0a6,l 

10 1 

2 3 ac (0b c ,o + 9 bc , e ) 



°i °oi °i 

i, r ii —0 0oi.i) r ia 



1 >fa 

2 (0Ol,a - 01a,( 



10 

--0 ac (0 



1 



01, c 



0ic,e), 



1 o 



aC 06c.li 





t-i a 
1 be 



) 1 01 

1 1 

2 ac 0fe c ,e, 

1 °ad/° N 1 *0a 

2 (06d,c + 0ca-,fe - 0f> c ,d) + - 0° a 06c,, 



The nonzero components of the Ricci curvature at the orders e 2 and e 1 



arc 



-2 1 ab 1 °ac °6d 

ii 00= ~2^ a 9 a b,e),e — ^ 9 ac 9bc,e9 9 a d,e 
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1 O m 0.0 

+ 2 9 %i,e9 ab 9 ab , e , (5.24) 

Roi= -(9 01 Soi,l)fl " \{9 ab 9 a b,i)e ~ ~ 9 ac 9 bc>6 9 bd L,i, (5-25) 
-l o 01 fa 1 0,0 

-Rafc = — 5 [9 a b,W — 7^ ^ \9acfi9bd,l + 9 a c,l9bdfi) 
10,0 \ 

+ 4 9 cd (9cd,l9 ab ,e + 9 cd ,e9ab,i) J ■ (5.26) 

-1 1 1 10 ac O l fe 10 01 1 °ab° 

i?oo= #a,06> ~ 2 ^ 9bc ' e9 a,e + 2 9 9oi y e9 a a fi ~ (9 a 9 a b,e)o 

1 0, , 1 O m , 

-- 9 ac 9 bc ,99 bd 9ad,o +- 9 m 9 ab (9oi,e9abfi + 9 i fi 9 ab ,e) (5.27) 

-1 1 O nl l n . 1 , 00 m 1 n , 

Roa= ^(9ab9 m (9oi9 Ub )e)e + j 9 cd 9 c d,e9 a b9 m (9 Q i9 m )e 

Moi °cd° ^ l , Q bc \ 1 ° bc de 

--(5 U1 9oi, a + 9 cd 9cd,a)e + t;(9 9 a bfi)c + j 9 bc 9 a bfi9 de 9 de , c 



1 O nl 0,0 lo,, 

+ 4 5 U1 9m A 9 cd 9cd,e 9 bd 9 c d,e9 ce 9 be ,a, (5.28) 



6 The asymptotic expansion 

We choose coordinates 

(x°, x 1 , x 2 , x 3 ) = (u, v, y, z) (6-1) 

in which the metric has the form ( |5.1| ). We introduce functions M, U, V, 
W of (9,v,y,z) such that 

5oi = -e~ M , 
o \ _ / e-^coshTy -e-^sinhP^ \ 
9ab ~ [ -e-^sinhiy e- u - v coshW ■ ^ 



It follows from (5^), (64), and (£^2) that the leading order metric has the 
form of the colliding plane wave metric, 

g= -2e~ M du dv + e~ v \e v cosh Wdy 2 - 2 sinh Wdy dz 

+e~ v cosh Wdz 2 ). (6.3) 

^From ( [5. 24 ), the only component of the leading order perturbation 
equation ( f4.13| ) which is not identically satisfied is 

Roo= 0. (6.4) 
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Using ( [5.24 ) and Q6.2| ) in ( |6.4| ), we obtain the ^-constraint equation, 
= I ( U + v e cosh 2 W + Wf) - UgMg. 



(6.5) 



^From (|5.25] )-( |5T2"8| ), the only components of the first order perturbation 
equation Q4.14D which are not identically satisfied are 



Roi= 0, Rab= 0, 
i?oo= 0, R 0a = 0. 



(6.6) 
(6.7) 



Using ( 5.25 )-( 5.26 ) and ( |6.2| ) in ( |6.6| ), we get the evolution equations in the 
colliding plane wave equations, 



Ug v 
V ev 

W 6v 



UgU v , 

i (U e V v + U v V e ) - (V e W v + V v Wg) tanh W, 
- {UgW v + U v W e ) + VgV v sinh W cosh W. 

i (-17*17,, + VgV v COsh 2 W + WflW^) . 



(6.8) 
(6.9) 

(6.10) 

(6.11) 



^From ( p^) and ( 5.27| )-( |5.28| ), we find that fl6.7|) is satisfied by a suitable 

l l 

choice of the first order metric components 9 a b, 9i a - 



7 Variational principle 

The variational principle for the vacuum Einstein field equations is 



5S = 0, 



S 



Ld x, 



L = Ry/-detg, (7.1) 
where R is the scalar curvature, 

R = g a ^R a p- 

Using ([4.2;), (|4.3|), and ( 4.10| ) to expand the scalar curvature, we obtain 



1-2 l-i 

R= — R +- R +0(1), 

e z e 

-2 a -2 

R=9 al3 R a/3 , 

-1 a -1 1 a -2 



(7.2) 
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For a metric of the form ( |5,21 ), we find that 



-2 

R=0, 

-10,-1 nn -1 l nn -2 

R=9 ab Rab +2 9 01 Roi - 9 00 R 00 . (7.3) 



The only order one metric component which appears in (7.3) is 

i , 



A = -g 



00 



In the derivation of the asymptotic equations, we used a coordinate system 
in which A = — see ( |5.3[ ). In the variational principle, A acts as a Lagrange 
multiplier for the constraint equation, so we do not set it to zero until after 
we take variations. 

We use ( |7.3| ) in (|7.1|), expand the result with respect to e, and write the 
expanded Lagrangian in terms of A and the functions M, U, V, W, defined 
in (|BT^). This gives 

L = - L +0(1), 



with 



L= {-2M 8v - 4U 9v + 3U e U v + V e V v cosh 2 W + W e W v } e~ u 



+A \u ee - ~ (U$ + V e 2 cosh 2 W + W$) + U e M e \ e 



M-U 



(7.4) 



We make the change of variables in the integration 

d 4 x = du dv dy dz = edO dv dy dz, 

and neglect the integration with respect to the parametric variables (y,z). 
The leading order asymptotic variational principle then becomes 

o o f -l 

8S=0, S= L dOdv. 



o 

Variations of S with respect to the first order metric component A give the 
constraint ( |6.5| ). Variations with respect to M, U, V, W give the evolution 
equations ( |6 . 8[ ) ( |67TT ) , after we set A = 0. It is permissible to set A = 



because the constraint is a gauge-type constraint which is preserved by the 
evolution equations. 
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8 Boundary conditions 



In this section, we discuss the derivation of boundary conditions for the 
asymptotic equations. For simplicity, we consider a "sandwich" wave located 
near the null surface u = which varies rapidly in a thin strip 

9- < - £ < 9+. 

We denote the slowly varying metrics on either side of the wave by 

J g+ in u > , . 

g " \ g in u < • (8 - 1} 

We consider a coordinate patch around a point on the surface u = with 
local coordinates (u, v, y, z) chosen as in the derivation of the asymptotic 
solution. In order for the metric outside the wave to join continuously with 
the solution inside, we must have 

g± -► -2e~ M± du dv + e~ u± (e v± cosh W±dy 2 

-2smhW±dydz + e~ v± cosh W±dz 2 ), (8.2) 



as u — > ± , where M±, U±, V±, W± are functions of (v,y,z). From (j6 
( S^2|), and the continuity of the metric, it follows that the solution of ( |6.8D - 
( 6,11| ) must satisfy the characteristic boundary conditions, 

M = M±, U = U±, V = V±, W = W±, when 6 = 6±. (8.3) 

This data need not satisfy the constraint ( |2.13| ). 

The asymptotic equations must be supplemented by a condition which 
specifies the profile of the wave. For example, we can impose a characteristic 
initial condition 

M = M , U = U , V = V , W = W , whent; = 0, (8.4) 

where Mq, Uq, Vq, Wq are functions of (9,y,z) which satisfy the constraint 
( |6.5D . The characteristic initial data must also be compatible with the char- 
acteristic boundary data, meaning that 

M (e ± ,y,z)=M±(p,y,z), 

together with analogous conditions for the other variables. 

Equations ( |6.8| ) — ( JBTTT ) , the characteristic initial condition ( |S.4| ) on v = 0, 
and the characteristic boundary condition (|8.3|) on 9 = 9_ form a well-posed 
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problem. Provided that the solution inside the wave is free of singularities, 
this problem has a unique solution. In particular, the solution at 9 = 9 + is 
uniquely determined. Thus, in principle, the asymptotic equations, together 
with the characteristic initial data ( |8.4|) , determine a set of jump relations 
which connect the minus and plus metrics ahead of and behind the wave, 
respectively. If the metric ahead of the wave is known, then the jump condi- 
tions provide characteristic boundary conditions on u = for the space-time 
behind the wave. Together with a characteristic initial condition on v = 
and u > 0, for example, this gives a characteristic initial value problem 
|0] for the full field equations. This problem determines the slowly varying 
metric behind the wave (at least locally). 

For instance, in the case of a plane polarized wave, the solution of ( |3.9| ) 
for U is g] 

U(9,v) = -log[f(9)+g(v)]. (8.5) 

Here / and g are functions of integration, and we do not explicitly show the 
possible parametric dependence of the functions on (y, z). From (|B.3|), (|87 
and (|8.5|) we have 

f(9) + g(0) = e~ u ^\ f(9-) + g(v) = e~ u -^. 

The solution is nonsingular provided that f(9) + g(v) > 0. 
It follows from (|8.5|) that the jump in U satisfies 



e -U+(v) _ e -V-(v) _ e -U (6+) _ e -t/o(0-). 
Use of ( |8.5| ) in ( |3.10| ) gives a linear wave equation for V, 
U + 9)V ev = \(g v Vg + f e V v ). 

Solution of this equation with the characteristic initial data V = Vq on 
v = and the characteristic boundary data V = V— on 9 = 9- determines, 
in principle, the solution V = V+ on 9 = 9 + . Finally, when W = 0, we 
define the f-constraint function G by 

G = U vv — ^ (U% + V%) + U V M V . (8.6) 

It follows from Q and (U)-(|l|) that 



Gq = UqG. 
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Integration of this equation with respect to 6 implies that 
logG + (v) -logG_(u) = U + {v) - U-(v). 

This equation provides a jump condition for M. 

One difficulty which arises in the formulation of boundary conditions 
ahead of the wave is that the metric g_ may not be given in a coordi- 
nate system which is compatible with the coordinate system used in the 
derivation of the asymptotic equations. It is then necessary to construct 
compatible coordinates (u,v,y,z). The u-coordinate is the phase, so it is 
a null coordinate of the metric which can be found by solving an eikonal 
equation, subject to appropriate initial conditions. The u-coordinate is a 
null coordinate which is orthogonal to u, while the y and z coordinates 
parametrize the null geodesies on the surface u = v = 0. 

If the gravitational wavefront u = forms a caustic, then the solution of 
the eikonal equation becomes multi-valued. When this happens, the local 
plane-wave approximation breaks down, and the asymptotic solution is not 
valid. However, the focusing at a caustic of the congruence of null geodesies 
associated with the phase does not necessarily imply the formation of a 
space-time singularity. 

9 Examples 

In this Section, we derive boundary conditions for the asymptotic equa- 
tions which correspond to the propagation of a non-planar gravitational 
wave into Minkowski space-time, the exterior Schwarzschild space-time, and 
Robertson- Walker space-time. In each example, we consider the case of 
spherical waves, where the boundary data can be explicitly computed. In 
this paper, we do not attempt to explore the physical consequences of the 
asymptotic equations in any detail. Our aim here is simply to illustrate how 
to apply the asymptotic equations to specific physical problems. 

9.1 Nonplanar wave propagation into Minkowski space-time 

We suppose that the space-time ahead of the wave is flat. In inertial coor- 
dinates (t,x), with t = x° and x = (a: 1 , a; 2 , a: 3 ), the metric is 

g = -dt 2 + dx 2 . 

We consider a wave with phase 

t — w(x) 

u = -= — . 

V2 
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The phase u is a null function of g_ if 

\Vw\ 2 = 1, 

where V is the gradient with respect to x. We define 

t + w(x) 

and choose coordinates y(x), z{x) such that Vro, Vy, Vz are orthogonal. In 
the (u, v, y, z) coordinates, we have 



A comparison of (|9.l| ) with (|8.2| ) shows that the minus boundary data is 
given by 



-u- _ 1 _ |Vy| 



M_ = 0, e - "- = — ,,_ ,, , e 



|Vy||Vz|| u=0 |Vz| 



, W- = 0. 

«=o 



For example, in the case of an outgoing spherical wave, suitable coordi- 
nates are 

t — r t + r 

u = —f=-, v = —j=-, y = V, z = f, (9.2) 

where (r, ip) are spherical polar coordinates and t > 0. In (u, t>, y, z) 
coordinates, the flat space-time metric is 

g = -2dudv + -(u— v) 2 (dy 2 + sin 2 y dz 2 ) . 

Evaluation of this metric at u = and a comparison with ( |8.2| ) gives the 
minus boundary data 

M_ = 0, e _c/ - = ^v 2 siny, e" y - = siny, W- = 0, 

where u > 0. In this case, M_, V_, and W- are independent of v, while 
U = U- satisfies the equation 

U - -U 2 

Thus, the boundary data satisfies the v -constraint equation ( |2.13| ). The 
solution is therefore identical to an exact solution for the collision of outgo- 
ing and incoming spherical waves, with the additional possibility of a slow 
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parametric dependence on the polar angles (•&, (p). Some exact solutions for 
spherical wave propagation into flat space-time are constructed in ||. 
For an incoming spherical wave, we use 

t + r t — r 

where t < 0. This leads to the same boundary data as in the case of an 
outgoing spherical wave, but with v < instead of v > 0. 



9.2 Gravitational waves incident on a black hole 

The exterior Schwarzschild metric is 

g = -adt 2 + -dr 2 + r 2 id® 2 + sin 2 $d(p 2 ) , (9.3) 
a 

where r > 2m and 

2m 

a{r) = 1 . 

r 

The contravariant metric tensor is 

gi = --d 2 + ad 2 + \(d 2 +'dl 
a r z \ sin v Y 

For simplicity, we consider an axially symmetric phase of the form 

t — w(r, 
u = -= . 

The function u is null if w satisfies the eikonal equation 

2 1 2 1 
aw r + -zw$ = -. 

r z a 

We define the orthogonal null coordinate 

t + w(r, 

v = ^^' 

and choose a coordinate y(r, ■d) whose gradient is orthogonal to the gradient 
of w(r, In that case, we have 

Vr = — ^2~' y$ = haw r, 
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where h(r, i?) is a suitable integrating factor. We take z = <p. In (u, v, y, z) 
coordinates, the Schwarzschild metric fl9,3|) is given by 



g = -2adu dv + -r^dy 2 + r 2 sin 2 fldz 2 . (9.4) 
h z 

Inversion of the change of coordinates (t, r, i9, ip) \— > (u, u, y, 21) implies that 
r = r-(v, y) and i? = ~d-(v, y) on u = for suitable functions r_ and A 
comparison of (9.4) with ( |8.2| ) implies that the boundary data is given by 



rr r 2 sin??- ,/ 

e" M -=a_, 6"^- = -— , e" y - = /i_sini?_, W_ = 0, 

rt_ 

where a_ = a(r_) and /i_ = /t(r_,$_). 

In the case of an incoming spherical wave incident on the black hole, 
suitable coordinates are 

t + A{r) t-A(r) 
u = — — ^, v = — -J^—i y = $i * = <P, (9-5) 

where 

a 

Integration of this equation implies that 

A(r) = r + log(r — 2m). 

In (u,v,y,z) coordinates, the exterior Schwarzschild metric is 

g = —2a dudv + r 2 (dy 2 + sin 2 y dz 2 ) . (9.6) 

^From Q9.5| ), we have r = r-(v) on u = where the function r_(v) is the 
solution of 

A(r_) = ~. (9.7) 



A comparison of ( |9.6| ) with ( |8.2j ) implies that the boundary data ahead of 
the incoming spherical wave is given by 



e~ M - = a-, e~ u - = r± siny, e~ v - = siny, W- = 0. (9.8) 

Dropping the minus subscripts, we find that the constraint function G in 
(|8.6|) for the boundary data fl9.g|) is given by 



q 2 ' ®"oTv r vv 

ar r 
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Differentation of (9.7) with respect to v implies that 

a a v 
r v = -- 7 =, r vv = - — . 

Use of this equation in the expression for G implies that G = 0. Thus, the 
boundary data ( |9.8|) satisfies the w-constraint equation ( 2.13| ), 



Numerical solutions of the interaction of a spherical gravitational wave 
with a black hole appear in ||]. 

9.3 Gravitational waves in a Robertson- Walker space-time 

The Robertson- Walker metric is 

g_ = -dt 2 + i { Y^dr 2 + r 2 (cM 2 + sin 2 #d<p 2 ) } , (9.9) 

where R(t) is the scale factor, and k = —1,0, 1. 

As in the Schwarzschild example, we consider an axially symmetric phase 
for simplicity, given by 

I(t)-w(r,$) 

where 

It = R, (1 - kr 2 ) w 2 + -^wl = 1. (9.10) 
We define an orthogonal null coordinate v by 

I(t) +w(r,&) 

v ~ 71 ' 

We choose a coordinate y(r, •&) whose gradient is orthogonal to the gradient 
of w(r, so that 

hw$ 



r 2 \/l — kr 2 ' 



y& = hy 1 — kr 2 w. 



where h(r, i?) is a suitable integrating factor, and take z = ip. In (u, v, y, z) 
coordinates, the Robertson- Walker metric (|9.9|) is given by 



2 V*^ 

g = --j^dudv + -jip^dy 2 + sin 2 -ddz 2 . (9.11) 

A comparison of (|9.11| ) with (|8.2| ) implies that the boundary data is given 
by 

e~ M - = ^T, e- u - = r ^^, e- v -=h„ S m^, W„ = 0, 
ri_ tl—TL_ 
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where r = r_(u, y), i? = $-(v, y), R = R-(v, y), and h = h-(v, y) on u = 0. 

For an outgoing spherical wave in a Robertson- Walker space-time, suit- 
able coordinates are 

I(t) + w(r) 



I(t) - w(r) 



where 



V2 
1 



y = #, z = tp, 



y/1 — kr 2 

Integration of this equation implies that 

sin -1 r if k = 1, 
w(r) = < r if k = 0, 

sinh -1 r if A; = — 

V 

The corresponding boundary data is given by 



-M- 



1 



rising ^_ v _ 



siny, W- = 0, 



(9.12) 



lit-) 



(9.13) 



" R 2 _ ' i? 2 
where t^(v) and r_(u) are given by 

sin(t;/V2) iffc = l, 
v/y/2 if jfe = 0, 

sinh(u/>/2) if fc = -1, 

and R_ = fl(t_). 

Dropping the minus subscripts, we find that the constraint function G 
in ( p^q ) for the boundary data ( |9.12j ) is given by 

7/ "" '' r,A (9.14) 



G 



i? 



iFrom flOD and ( pTlSD , we find that 



2 ' ™~ 2i? 3 ' 



Use of these expressions in fl9.14Q gives 

G= RRu_Rl + k 

Thus, in general, the boundary data ( |9.12| ) does not satisfy the u-constraint 
equation (2.13). 
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